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§ : ABSTRACT 
(N 

H ' We present an analysis of the compact X-ray source population in the Seyfert 2 

galaxy NGC 1068, imaged with a ~ 50 ks Chandra observation. We find a total of 84 
compact sources on the S3 chip, of which 66 are located within the 25.0 B-magnitude 
(arc sec)^^ isophote of the galactic disk of NGC 1068. Spectra have been obtained 
for the 21 sources with at least 50 counts, and modeled with both multi-color disk 
blackbody and power-law models. The power-law model provides the better description 
\ of the spectrum for 18 of these sources. For fainter sources, the spectral index has been 

I estimated from the hardness ratio. Five sources have 0.4-8 keV intrinsic luminosities 

' greater than 10^^ erg s^^, assuming that their emission is isotropic and that they are 

associated with NGC 1068. We refer to these sources as Intermediate Luminosity X- 
^) • ray Objects (IXOs). If these five sources are X-ray binaries accreting with luminosities 

^ . that are both sub-Eddington and isotropic, then the implied source masses are ^ 7 Mq, 

and so they are inferred to be black holes. Most of the spectrally modeled sources 



> 



O I have spectral shapes similar to Galactic black hole candidates. However, the brightest 

I compact source in NGC 1068 has a spectrum which is much harder than that found in 

' Galactic black hole candidates and other IXOs. The brightest source also shows large 

amplitude variability on both short-term and long-term timescales, with the count rate 
^ ' possibly decreasing by a factor of two in ~ 2 ks during our Chandra observation, and 

^ . the source flux decreasing by a factor of five between our observation and the grating 

observations taken just over nine months later. The ratio of the number of sources with 
luminosities greater than 2.1 x 10^*^ erg s~^ in the 0.4-8 keV band to the rate of massive 
(> 5Mq) star formation is the same, to within a factor of two, for NGC 1068, the 
Antennae, NGC 5194 (the main galaxy in M51), and the Circinus galaxy. This suggests 
that the rate of production of X-ray binaries per massive star is approximately the same 
for galaxies with currently active star formation, including "starbursts" . 

Subject headings: accretion, accretion disks — galaxies: Seyfert — galaxies: individual 
(NGC 1068) — galaxies: starburst — X-rays: binaries — X-rays: individual (NGC 
1068) 
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1. INTRODUCTION 

Observations with the Einstein Observatory were the first to reveal the existence in some 
nearby galaxies of non-nuclear compact sources with X-ray luminosities^ exceeding 10^^ erg s~^, 
which is well above the Eddington limit for a 1.4 solar mass (Mq) neutron star (for a review of 
the Einstein observations see e.g., Fabbiano 1995). More recent observations with ROSAT have 
shown that these Intermediate Luminosity X-ray Objects (IXOs)'^ are common, with 9 of the 29 
galaxies surveyed by Lira, Lawrence, & Johnson (2000) hosting at least one source (see also Colbert 
& Mushotzky 1999 and Roberts & Warwick 2000 for comparable studies with ROSAT). 

If these sources are powered by accretion, and are emitting at close to their Eddington limit, 
then they would contain black holes with masses up to several hundred Mq or higher (e.g., Zezas, 
Georgantopoulos, &; Ward 1999; Kaaret et al. 2001; Matsumoto et al. 2001). Black holes of this 
mass cannot arise from the collapse of a single, non-zero metallicity star, but may have instead 
originated from an earlier generation of zero metallicity (population III) stars (e.g., Madau &: Rees 
2001). Alternatively, a massive black hole may have formed in the center of a globular cluster. If 
the cluster then merges with the galaxy disk, the black hole may then accrete gas in dense molecular 
clouds (Miller & Hamilton 2002). Ebisuzaki et al. (2001) and Portegies Zwart k McMillan (2002) 
have suggested black hole formation in young star clusters by rapid collapse. 

If the X-ray emission is beamed, then the required black hole mass may be reduced to val- 
ues found in Galactic black hole candidates. Reynolds et al. (1997) and Kording et al. (2002) 
have suggested that some of the X-ray emission from IXOs may be relativistically beamed, and 
that IXOs may be objects similar to the Galactic superluminal sources such as GRS 1915-1-105 
and GRO J1655— 40, but with their jets pointed towards the observer. Alternatively, IXOs may 
represent a brief evolutionary phase of close separation X-ray binaries, in which the neutron star 
or black hole accretes mass from an early type star at a highly super-Eddington rate without the 
system evolving into a common envelope (King et al. 2001). This model invokes a thick accretion 
disk with a central funnel, which may radiate a total luminosity in excess of the Eddington limit 
(e.g., Abramowicz et al. 1980), although the stability of such structures is questionable. 

Other possibilities include supernovae exploding into dense circumstellar material (e.g., Franco 
et al. 1993; Plewa 1995). Such objects can reach luminosities of up to 10^^ erg (e.g., SN 1988Z: 

Fabian & Terlevich 1996; SN 1995N: Fox et al. 2000), similar to that observed in the most luminous 
IXOs. However, supernovae cannot show the kind of X-ray variability noticed in some IXOs (sec 
below). More recently, Begelman (2002) has argued that an accretion disk around a solar mass 
black hole could exhibit strong density inhomogeneities, which would allow the escaping flux to 



^Unless noted otherwise, all luminosities quoted in this paper assume isotropic emission. 

^Here and elsewhere in this paper, we refer to this class of object as Intermediate Luminosity X-ray Objects. 
Other acronyms for this class include Ultra-Luminous X-ray sources (ULXs), Super-Luminous X-ray Sources (SLS), 
and Super-Eddington X-ray Sources (SES). 
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exceed the Eddington limit by a factor of 10-100, thus explaining the high luminosities observed 
without beaming. 

Some of the recent work on the X-ray spectra of IXOs has concentrated on the contribution 
from an optically thick accretion disk around a black hole (e.g., Makishima et al. 2000). Although 
many of the IXO spectra can be modeled by a multi-color disk blackbody (MCD) continuum, the 
highest (i.e., innermost) disk temperature derived from the observations is usually larger than that 
expected for an accretion disk around a black hole accreting at, or below the Eddington limit. 
One possible resolution of this problem is that the black hole is spinning in the same direction as 
the disk, and hence the inner radius is closer to the black hole, and the temperature higher, than 
for a non-rotating black hole of the same mass (e.g., Zhang, Cui, &; Chen 1997; Makishima et al. 
2000). Alternatively, the high temperatures may indicate the presence of a slim (i.e., advection 
dominated, but optically thick) disk, in which substantial X-radiation is emitted from within the 
last stable circular orbit around a Schwarzschild black hole (Watarai et al. 2000; Watarai, Mizuno, 
&: Mineshige 2001; Mizuno, Kubota, & Makishima 2001). 

More recent observations with Chandra have identified IXOs that cannot be modeled with an 
MCD spectrum, but require instead a hard power-law continuum of photon index T 1-2 (e.g., 
Kubota et al. 2001; Strickland et al. 2001). While some models of supernovae exploding into a 
dense circumstellar medium do predict hard power-law continua with T ~ 1.6-2 (Franco et al. 1993; 
Plewa 1995), there is evidence of a spectral transition between a high (with a soft X-ray spectrum) 
and a low (with a hard X-ray spectrum) intensity state in two IXOs in IC 342 (Kubota et al. 2001), 
and in one IXO in the dwarf galaxy Holmberg IX (La Parola et al. 2001). These observations lend 
further weight to the X-ray binary interpretation of IXOs, in which the hard power-law continuum 
is interpreted as evidence for either a low (hard) state, as commonly seen in Galactic black hole 
candidates, or highly Comptonized, optically thick disk emission, which often characterizes the very 
high state of Galactic black hole candidates (e.g., Kubota, Done, &: Makishima 2003). Conversely, 
one source in M51 exhibited hard (F = 1.2), luminous (L(0.5 — 8keV) = 2.7 x 10^^ erg s~^) emission 
in June 2000, but one year later showed a very soft (F > 5.1), lower luminosity (L(0.5 — 8keV) = 
5.6 X 10^^ erg s~^) spectrum (Terashima & Wilson 2002a). Such behaviour is reminiscent of soft 
X-ray transients. 

Variations in the X-ray emission of several IXOs have been observed both during a single 
observation (e.g., Okada et al. 1998; Zezas et al. 1999), and between observations which span 
several years (e.g.. La Parola et al. 2001). However, the most spectacular evidence for X-ray 
periodicity in an IXO has been found by Bauer et al. (2001) for a bright compact source in the 
Circinus galaxy (CG X-1 or CXOU J141312.3-652013; see also Smith & Wilson 2001). The period 
is 27.0 ± 0.7 ks, and Bauer et al. (2001) argue against the notion that this source is a foreground 
Galactic AM Her type system based primarily on the low surface density of Galactic X-ray sources 
in this direction (/ = 311°, 6 = — 3.°8). Another IXO in IC 342 is possibly periodic, with a period 
of 31 or 41 hours, in agreement with that expected for a semi-detached binary consisting of a black 
hole and a main-sequence star of tens of solar masses (Sugiho et al. 2001). 
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Constraints on the evolution and nature of IXOs can be obtained by comparing the X-ray 
luminosity functions (XLFs) of the X-ray source populations in various types of galaxies and in 
various evolutionary stages. For example, the slopes of XLFs for starburst galaxies (e.g, M82, "the 
Antennae") are flatter than those found in early-type galaxies (e.g., NGC 1553, NGC 4967), but 
are similar to that of the high-mass X-ray binary population in our own Galaxy (Kilgard et al. 
2002; Grimm, Gilfanov, & Sunyaev 2002; Zezas k. Fabbiano 2002). IXOs tend to be associated with 
regions of active star formation in spirals, but there are also IXOs in elliptical galaxies (Colbert k. 
Ptak 2002). This strong association with star formation suggests that the brightest X-ray sources 
in spiral and starburst galaxies are likely to be young, short-lived sources, e.g.. X-ray binaries with 
O and B type companions, or supernova remnants (Kilgard et al. 2002). In early- type galaxies 
(ellipticals and lenticulars), the XLF tends to break at the Eddington luminosity of a 1.4 M0 
neutron star (e.g., Sarazin, Irwin, &; Bregman 2000; Blanton, Sarazin, & Irwin 2001), but this 
seems not to be the case in active star-forming regions. 

The galaxy NGC 1068 is not only the most luminous (~ 10-^^ Lq) nearby Seyfert 2, but also 
one of the most luminous (also ~ 10^^ Lq) starbursts in the local universe (e.g., Telesco & Decher 
1988). It has long been suggested (e.g., Weedman 1983) that the two phenomena are related, but 
the exact process remains elusive. The starburst is circumnuclear, in the galaxy disk and on a scale 
of ~ 2 kpc. It is thus a very different phenomenon to many other regions of star formation in the 
universe, such as "grand-design" late-type spirals (e.g., M51) and star formation induced by galaxy 
mergers (e.g., the Antennae). For this reason, we felt it worthwhile to investigate the compact 
X-ray source population in a galaxy with such a luminous, circumnuclear starburst. In this paper, 
we adopt a distance of 14.4 Mpc to NGC 1068, so 1" = 70 pc (e.g., Bland-Hawthorn et al. 1997). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations and Initial Analysis 

NGC 1068 has been observed on four occasions with the Chandra X-ray Observatory (Weis- 
skopf et al. 2001) in direct imaging mode with the Advanced CCD Imaging Spectrometer (ACIS; 
Garmire 1997) at the focal plane of the High- Resolution Mirror Assembly (van Speybroeck 1997). 
Results from an analysis of the nuclear and extended X-ray emission associated with NGC 1068 
are presented elsewhere (Young, Wilson, &: Shopbell 2001, hereafter Paper I). A full account of the 
observations is given in Paper I. 

We concern ourselves here with a study of the discrete X-ray source population in NGC 1068, 
imaged within the 8.'4 x 8.'4 (35.3 x 35.3 kpc) field of view of the S3 chip. The data reduction and 
analysis were done using the Chandra Interactive Analysis of Observations (CIAO) software version 
2.2.1 (released on December 13, 2001) and CALDB 2.12 (released on February 14, 2002), and the 
reprocessed (on March 10, 2001 using CALDB 2.3 and version R4CU5UPD14.4 of the processing 
software) event files. New level 2 event files were created, applying the latest telescope geometry 
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and detector gain corrections, and including the same ASCA grades, bit status, and time filters as 
in the existing level 2 event files. Periods of high and low background (i.e., flares or data dropouts 
due to telemetry saturation) were excluded from the data. This was achieved by creating a light 
curve over the full energy range for the whole S3 chip, excluding the brightest sources of X-ray 
emission, and removing events ±3(t from the mean count rate. This procedure gives an effective 
exposure time of 47.1 ks (corrected for the dead-time in the detector). We have checked whether the 
decline in low energy quantum efficiency of the S3 chip affects our results by using the ACISABS 
spectral model on a few sources, but found the spectra to be insignificantly different, presumably 
because our observations were taken early in the mission (2000 Feb. 21). 

2.2. Source Detection 

An image from the 3.2s frame time observation of NGC 1068 (obsid 344), on the scale of 
the galactic disk, is shown in Fig. 1. There is a considerable amount of diffuse X-ray emission in 
this image, which extends at least 60" (4.2 kpc) to the northeast, 50" (3.5 kpc) to the southwest, 
20" (1.4 kpc) to the northwest, and 30" (2.1 kpc) to the southeast of the nucleus. The dominant 
large-scale structures arc "spiral arms" that curve to a lower position angle (P.A.) with increasing 
galactocentric distance (Young et al. 2001). Many compact sources of X-ray emission associated 
with NGC 1068 are also seen in this image. The luminous X-ray sources appear to be concentrated 
in the starburst region, but there are also many sources located further out from the center of the 
galaxy. 

Images were extracted from the reprocessed level 2 events file in soft (0.4-1.5 keV), hard (1.5- 
5.0 keV), and full (0.4-5.0 keV) energy bands. We have used the CIAO program wavdetect to search 
the images in the three energy bands for discrete sources of X-ray emission. This program uses a 
scalable Marr, or "Mexican Hat," wavelet function to parameterize the shape and extent of each 
source, and generally performs better in crowded fields than the standard "sliding cell" algorithm 
which has traditionally been used to detect sources in X-ray images (Freeman et al. 2002). We have 
analyzed the images using wavelet scales in the range 1 pixel (0."492) to 16 pixels (7."87), separated 
by a factor of \/2 (i.e., 1, \/2, 2, 2\/2, 4, 4\/2, 8, 8\/2, and 16 pixels). These scales correspond to 
the radius of the Mexican Hat function, and are representative of the range of source sizes seen 
in each of the images. The wavelet source detection threshold was set to 10~^, which will give 
approximately one false source for the whole S3 chip.^ The total number of sources detected by 
wavdetect in the soft, hard, and full energy band images was 115, 67, and 138, respectively. Each 
of these sources was examined carefully by eye, and only those 84 sources which appear compact 
to the eye are included in the source list (Table 1). 



*See §11.2 of the CIAO Detect Manual publicly available at 
littp://asc.liarvard.edu/ciao/download/doc/detect_html_manual/Manual.litml. 
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2.3. Spectra 

Spectra were extracted for those 21 sources listed in Table 1 with at least 50 counts (after 
background subtraction). We have chosen to ignore 4 sources close to, or in the "prong" of bright, 
diffuse X-ray emission extending in P.A. ~ 30° near the Seyfert nucleus, since their spectra are 
highly uncertain, due to contamination by the diffuse emission, and their signal-to-noise ratio is 
< 7. The shape and size of the extraction region was chosen to include as much of the source 
flux as possible, while minimizing the contribution from background and nearby sources. Circular 
regions with diameters in the range l."8 to 6."0 sufficed for sources within ~ 2' of the telescope 
mirror axis. At off-axis angles greater than this, there is considerable broadening and distortion of 
the point spread function and so elliptical regions with major axis diameters in the range 3."7 to 
16."6 and major to minor axis ratios typically 1.2—1.6 were used. Background spectra were usually 
accumulated from circular or elliptical annuli surrounding the source. In the case of sources close to 
the nucleus, where the diffuse extended emission is strongest, background spectra were occasionally 
taken from circular regions adjacent to the source. Sources which lie within any given background 
region were excluded when calculating the background. We created a response matrix for each 
spectrum using the CIAO programs mkrmf and mkarf and the latest calibration files (version 
N0002, released in August 2001). Response matrices are available for 32 x 32 pixel regions on the 
S3 chip. Differences between response matrices for adjacent regions are very small, and so, for 
any given source, we used the response matrix for the region of the chip where most of the source 
photons are located. In the case of sources which straddle a node boundary, we extracted a response 
matrix from each node, and then created an average response matrix using weights proportional to 
the background subtracted counts from each node. We have ignored the effects due to the secular 
decline in the ACIS quantum efficiency since the observations of NGC 1068 took place early in the 
mission, and the effects on the derived values of the column density are negligible (less than 10%) 
compared to the statistical error. Prior to performing the spectral analysis with XSPEC version 
11.1.0 (Arnaud 1996), we rebinned the source spectra so that there were at least 15 counts per bin 
(for sources with > 100 counts) or at least 10 counts per bin (for sources with < 100 counts), thus 
allowing use of the statistic. The results of the spectral analysis are given in Table 2. 



2.4. Timing Analysis 

Light curves were extracted in the 0.4-5 keV band for each of the 21 non-nuclear sources with 
sufficient counts for spectral analysis (see Table 2). The low count rates for most of the sources 
preclude an analysis of the light curves based on the statistic, which requires the data to be 
binned. To assess the probability of variability, we have, instead, compared the distribution of 
event arrival times for each source with that of a large source-free region of the S3 chip (10557 
events), using the Kolmogorov-Smirnov (K-S) test (see e.g., Press et al. 1992, p. 623). 
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3. RESULTS 

3.1. Detected Sources 

Wc find a total of 84 compact sources on cliip S3 (Table 1), of which 66 are located within 
the 25.0 B-magnitude (arc sec)~^ isophote of the galactic disk of NGC 1068 (de Vaucouleurs et al. 
1991). Fig. 2 shows the locations of the 84 compact sources superposed on the Digitized Sky Survey 
image. As can be seen in this figure, approximately half of the sources are projected onto the bright 
inner disk region of NGC 1068. In order to quantify the expected number of background sources in 
the field, we have performed an identical analysis on an ACIS-S observation of 3C 47 (obsid 2129), 
which has a similar exposure time (41.6 ks) to our observation of NGC 1068. We have considered 
only regions > I' from 3C 47 as there is bright extended X-ray emission surrounding this source. 
Regions along the readout trail were also excluded from our analysis. We find a mean surface 
density of 0.37 sources (arc min)^^ in this observation, which is almost equal to the mean surface 
density of sources 0.45 sources (arc min)~^] outside the 25.0 B-magnitude (arc sec)~^ isophote 
of the galactic disk of NGC 1068. This suggests that most of the latter sources are background 
objects unrelated to NGC 1068. From the mean surface density of sources in the 3C 47 field, we 
would expect, in our Chandra observation, approximately 12 background sources projected onto 
the galactic disk of NGC 1068, whereas we detect 66 sources. This number of background sources 
is a conservative estimate, since the emission from such sources must be partially absorbed by gas 
in the galactic disk of NGC 1068. 

The expected number of background sources can also be estimated from the density of sources 
detected in the deep imaging sky surveys with Chandra (see e.g., Cowie et al. 2002 and references 
therein). In order to determine the limiting flux level during our Chandra observation, we need to 
compute the correlation of the Mexican Hat function with a two-dimensional Gaussian (representing 
the point spread function) using the procedure described in Appendix A of Freeman et al. (2002). 
The result of this correlation is then compared with the value expected for a given significance 
and background level (see Appendix B of Freeman et al. 2002). For sources beyond the 25.0 B- 
magnitude (arc sec)~^ isophote of the galactic disk, the point spread function is large, and we adopt 
a Gaussian of width a = 3.0 pixels, which is appropriate for a source 5' off-axis and an energy of 
4.5 keV. We find the detection threshold to be ~ 13.5 counts for our significance and background 
level. This number of counts corresponds to a unabsorbed flux of 1.8 x 10~^^ erg cm~^ s~^ in the 
0.4-5 keV band or 2.1 x 10~^^ erg cm~^ s~^ in the 2-8 keV band, assuming a power-law of photon 
index F = 1.2 and the Galactic column density towards NGC 1068 of iVH(Gal) = 3.53 x 10^° cm-^ 
(Dickey &; Lockman 1990). At this flux level, we would expect a surface density of 0.38 sources 
(arc min)~^ from the results of Chandra blank field observations in the 2-8 keV band (Cowie et al. 
2002). While this number is similar to the surface density of sources outside the 25.0 B-magnitude 
(arc sec)~^ isophote in the NGC 1068 observation, the comparison should be considered as only 
approximate, since the energy band (2-8 keV) used by Cowie et al. (2002) is different from the one 
used in our analysis (0.4-5 keV). 
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3.2. Spectral Analysis 

The source spectra were initially modeled in the 0.4-8 keV band with two alternative descrip- 
tions of the continuum: (i) a multi-color disk (MCD) blackbody spectrum,^ which represents the 
emission expected from an optically thick accretion disk (Mitsuda et al. 1984; Makishima et al. 
1986) and (ii) a power-law spectrum (photon spectral index T). We have included in the model 
absorption by both the Galactic column density towards NGC 1068 of A'^H(Gal) = 3.53 x 10^*^ cm~^ 
and a column density N}i intrinsic to the source; the atomic cross-sections and abundances for 
the absorption columns were taken from Morrison &: McCammon (1983) and Anders & Grevesse 
(1989), respectively. The results are given in Table 2. The MCD and power-law models both pro- 
vide a reasonable description of the data for each of the sources, with neither model being rejected 
at > 99% confidence using a test. Most of the source spectra have color temperatures and inner 
disk radii in the range fcTin ~ 0.5-2 keV and Ri^ ~ (10-200) x (cos^)"-*^/^ km, respectively, and 
intrinsic column densities up to A'^h ~ 6.0 x 10^^ cm~^ (MCD models), or photon spectral indices in 
the range F ~ 1-3 and intrinsic column densities up to Nn « 9.0 x 10^^ cm~^ (power-law models). 
We note that the best-fit intrinsic column densities tend to be larger for the power-law model than 
for the MCD model. It is also notable that 18 of the 21 sources in Table 2 are better modeled (i.e., 
have a smaller x^/^j where u is the number of degrees of freedom) by a power-law than an MCD 
model. Moreover, when the total set of spectra is considered, the MCD model can be rejected at 

> 99% confidence using a test. 

The two most luminous objects have contrasting spectral shapes (Fig. 3): the brightest 
source (J024238. 9— 000055) has a very hard spectrum (F = 0.85+°:n) 

whereas the other source 

(J024237.9— 000118) has a much softer spectrum (F = i.J^Qg). A significant improvement (at 

> 90% confidence using the F test for two additional free parameters) in the fit to the spectrum of 
J024240. 7— 000055 is obtained when we include a narrow Gaussian emission line. The line energy 
and normahzation are 1.031]':o;q3o keV and (1.02+Q;g^) x 10^^ photons cm~^ s~\ respectively. The 
equivalent width is 280^^gQ eV with respect to the observed (absorbed) power-law continuum, and 
the emission line energy is consistent with that of Ne X Lya. Finally, there is an absorption feature 
at ^ 1.6 keV in the spectrum of J024238. 9— 000055 (see Fig. 3), which may result from an artifact 
in the current release of the gain tables.^ 

In order to study the population statistically, we have adopted the following approach. For 
the sources whose spectra we have analyzed, the luminosities in the 0.4-8 keV band were estimated 
from the results of fitting a power-law continuum model (see Table 2). For the remaining sources. 



■''The MCD model in XSPEC returns the color temperature of the inner accretion disk kT^ and an inner disk 
radius rin = K^^^ D {cos9)~^^^ , where K is the normalization of the MCD spectrum, D is the distance to the source, 
and 9 is the (unknown) inclination of the disk (face on corresponds to ^ = 0°); the luminosity emitted from the disk 
into our line of sight varies as rf^^cosO. The true inner disk radius is _Rin = ^K^nn, whore ^ = 0.412 (Kubota et al. 
1998) and k — 1.7 (e.g., Shimura & Takahara 1995) is the ratio of the color temperature to the effective temperature. 

®http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_projects/index.html 



-9- 



which were spectrally unmodeled, we have assumed a power-law continuum model modified by (i) 
the Galactic column density plus (ii) an intrinsic column density of iVn = 1.28 x 10^^ cm~^. The 
latter column density is the average (e.g., Bevington & Robinson 1992, p. 58) of the spectrally 
measured intrinsic column densities (power-law model; Table 2) weighted by l/o"^, where a is the 
larger of the positive and negative errors on the intrinsic column density. We believe that these two 
values of column density are broadly representative of the range of column densities through which 
the X-ray sources in NGC 1068 are seen. From the hardness ratio HR = C(l. 5-5.0 keV)/C(0.4- 
1.5 kcV), where C represents the number of counts in the given band, we then determined the 
intrinsic photon index of each spectrally unmodeled source. These photon indices were then used 
to estimate the luminosities in the 0.4-8 keV band. In cases in which only an upper (lower) limit to 
the hardness ratio was determined from the data, we used this upper (lower) limit to determine the 
photon index. Figs. 4a and b are plots of HR against absorption-corrected luminosity. For sources 
too weak for spectral analysis, the Galactic column density is used in Fig. 4a and the Galactic plus 
average intrinsic column is used in Fig. 4b. It is notable that essentially all sources with 2 X 10^8 
erg s~^ < L(0.4 - 8keV) < 5 x 10^^ erg s"^ have hard spectra, with -0.5 < F < 2.5 (Fig. 4b). 

Under the assumption of absorption by the Galactic column density only, 5 and 15 sources 
have luminosities (after correcting for absorption) exceeding 10^^ and 5 x 10^^ erg s^^ in the 0.4- 
8 keV band, respectively. Alternatively, if an intrinsic column density of A^h = 1.28 x 10^^ cm~^ is 
assumed, then 5 and 18 sources have luminosities (after correcting for absorption) exceeding 10^^ 
and 5 x 10^^ erg s~^ in the 0.4-8 keV band, respectively. 

3.3. The Distribution of Source Luminosities 

The distribution of luminosities of the sources that we detect on chip S3 in the NGC 1068 
observation is shown in Fig. 5. We have excluded from this figure the three sources that were 
undetected in the 0.4- 5 keV band image (marked by open circles with upward and rightward 
pointing arrows in Fig. 4). Also shown in Fig. 5 is the distribution of luminosities (assuming their 
distances to be that of NGC 1068) for the sources that we detect on chip S3 in the 3C 47 observation 
(see § 3.1); we consider these sources as representative of the unrelated foreground and background 
objects in the NGC 1068 field. Most of the sources in the 3C 47 field have insufficient counts for a 
spectral analysis, and so we have assumed a power-law continuum (of photon index F = 2.0, typical 
of an AGN) modified by the Galactic column density towards 3C 47 of Ah (Gal) = 5.87 x 10^° cm~^ 
(Dickey & Lockman 1990). This spectrum was then used to estimate fluxes in the 0.4-8 keV band 
from the number of counts observed in the 0.4-5 keV band. We detect two sources in the soft 
(0.4-1.5 keV) band image of 3C 47 that were undetected in both the full (0.4-5 keV) and hard 
(1.5-5 keV) band images, and so we have used the number of counts observed in the soft band to 
estimate the flux in the 0.4-8 keV band. For luminosities greater than ~ 10^^ erg s~^, the surface 
density of sources in the 3C 47 field is much lower than that towards NGC 1068, so we can ignore 
the contribution of unrelated foreground and background objects to the luminous source population 
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in NGC 1068. There is, however, one source in the 3C 47 field with a luminosity of > 10 erg s~ 
in the 0.4-8 keV band at the distance of NGC 1068. This source has a flux of 1.3 x 10~^^ erg 
cm^^ s~^ in the 2-10 keV band (assuming a power-law of photon index T = 2.0 and the Galactic 
absorption towards 3C 47 of A^H(Gal) = 5.87 x lO^'' cm~^). At this flux level, we would expect 
~ 0.2 sources in the 3C 47 field, from the results of the ASCA Medium-Sensitivity Survey (Ueda 
et al. 1999), and so there is a ~ 20% chance of detecting a foreground or background source this 
bright in the NGC 1068 field. 

The observed distribution of source luminosities in NGC 1068 is subject to uncertainty due to 
the Poissonian nature of observed source counts, source confusion, incompleteness, and contamina- 
tion by sources that are spurious. While Monte Carlo simulations can be used to estimate the effect 
of these biases (e.g., Zezas &: Fabbiano 2002; Kim & Fabbiano 2002), this approach is unlikely to 
benefit us here, since there are 54 sources in Table 1 with a signal-to-noise ratio < 5 as a result of 
diffuse emission at small galactocentric distances in NGC 1068. We do not consider it worthwhile 
computing the cumulative XLF, since the error in the slope is likely to be much larger than the 
distribution of slopes presented in Kilgard et al. (2002). 

3.4. Source Variability 

The brightest source, J024238. 9— 000055, is possibly variable, with a 0.36% probability that 
the arrival times of the source events are drawn from the same parent population as those of the 
background events (and that the K-S statistic exceeds the observed value by chance). The light 
curve for this source is shown in Fig. 6. Although we cannot reject at > 99% confidence, using a 
test, the null hypothesis that the X-ray flux is constant, there is a suggestion of a rapid decrease in 
the source flux some 30 ks into the observation, where the count rate decreases by a factor of two 
in ~ 2 ks. 

In order to search for long term variability we compared our data with the Chandra grat- 
ing observations. Due to the much poorer signal-to-noise ratio in the grating data, we con- 
centrated on the two brightest sources. In the High Energy Transmission Grating Spectrome- 
ter (HETGS) observation performed on December 4-5, 2000 (obsid 332) for 46 ks, we detected 
49.9 ± 8.5 and 43.0 ± 7.7 counts (after background subtraction) in the 0.4-5 keV band from 
J024238.9-000055 and J024237.9-000118, respectively. We have also modeled the spectra of 
both sources with the corresponding power-law model given in Table 2, keeping all the param- 
eters fixed at their best-fit values except for the power-law normalization which was allowed to 
vary. The observed fluxes in the 0.4-8 keV band are 1.0 x 10^^'^ and 3.9 x 10~^^ erg cm~^ for 
J024238.9-000055 and J024237.9-000118, respectively These results suggest that the observed 
flux from J024238. 9— 000055 has declined significantly (by a factor of ~ 5) between the ACIS-S 
and HETGS observations, while the flux observed from J024237.9— 000118 has remained essentially 
unchanged. 
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The Low Energy Transmission Grating Spectrometer (LETGS) observation was performed 
~ 1 day later on December 5-6, 2000 (obsid 329) for 77 ks. Source J024237.9-000118 was outside 
the field of view of the zeroth order image, so we have only analyzed source J024238. 9— 000055. We 
detected 49.6 ±8.7 counts (after background subtraction) in the 0.4-5 keV band. A similar spectral 
modeling approach was adopted for the LETGS spectrum as that used for the spectrum from the 
HETGS observation. We find the observed fiux in the 0.4-8 keV band to be 1.1 x 10~^^ erg cm~^ s~^, 
similar to that observed in the HETGS observation. 

4. DISCUSSION 

4.1. The Intermediate Luminosity X-ray Objects (IXOs) 

We detect a total of 84 compact, non-nuclear X-ray sources on the S3 chip in our observation 
of NGC 1068. After correcting for only Galactic absorption, 5 of these sources have 0.4-8 keV 
luminosities greater than 10^^ erg s~^, assuming that their emission is isotropic and that they are 
located within NGC 1068. All of these sources project inside the 25.0 B-magnitude (arc sec)~^ 
isophotc of the galactic disk. In Fig. 7 wc show an optical continuum image of NGC 1068 at 6100 A 
(Poggc & dc Robcrtis 1993) with the positions of 4 of these sources superposed [the other source 
is located outside the 2'.0 x 2'.0 (8.4 x 8.4 kpc) field of view]. Also shown in this image are the 
positions of the sources which have luminosities (corrected for absorption) between 5 x 10^^ erg s"^ 
and 10^^ erg s~^ in the 0.4-8 keV band. There are 3 sources in this luminosity range outside 
the field of Fig. 7. Two of these project outside the 25.0 B-magnitude (arc sec)~^ isophote of the 
galactic disk.^ Their hardness ratios are HR = 2.0 ib 0.9 and HR = 0.33 it 0.06, respectively. For 
a source with a power-law spectrum of photon index F = 2.0, typical of an AGN, HR = 0.34 if we 
assume absorption from the Galactic column density only, and HR = 0.51 if a column density of 
Nu = 1.28 X 10^^ cm~^ is assumed. Thus, we conclude that the spectrum of one of these sources 
is consistent with that of a typical AGN, while the spectrum of the other source is consistent with 
that of a typical AGN only if it is highly absorbed. 

If the 5 sources with 0.4-8 keV luminosities greater than 10^^ erg are X-ray binaries 
accreting with luminosities that arc sub-Eddington (i.e., ^ 1.5 x 10^^ (M/Mq) erg s~^, where 
M and Mq are the masses of the X-ray source and the Sun, respectively), then the implied source 
masses are ^ 7-100 M©. Such objects are almost certainly black holes, since the upper limit on the 
mass of a non-rotating neutron star is ~ 3Mq (for a review see e.g., van Paradijs &: McClintock 
1995). The maximum mass of a black hole formed by the collapse of a single, non-zero metallicity 
star within the mass range ~ 25-40 M© is commonly believed to be ~ 15 Mq due to the effects of 



^Thoro arc 11 sources marked in Fig. 7 with luminosities between 5 x 10^* erg s~^ and 10^® erg s~^ (assuming 
an association with NGC 1068). We detect 10 sources in this luminosity range on chip S3 assuming only Galactic 
absorption, and 13 sources if both Galactic and internal absorption is assumed. 
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mass loss through stellar winds. The black hole mass can be even lower 10 Mq) for black hole 

progenitors in close separation binary systems since, prior to collapse, the progenitor's helium core 
may undergo mass loss in a Wolf-Rayet phase if the hydrogen layers are removed by the companion 
star in a common envelope phase (Fryer & Kalogera 2001). However, the limit on the black hole 
mass may be higher, since stars more massive than ~ 40 Mq may not produce a supernova explosion 
and all of the stellar material falls back onto the black hole (Fryer 1999). Moreover, the mass loss 
rates for low metallicity stars may be small, thus allowing the black hole progenitors to retain 
most of their original stellar mass (Vink, de Koter, &: Lamers 2001). Stars with masses within 
the range ~ 100-250 Mq may undergo explosive nuclear burning at the end of their evolution, 
however, and no remnant is left (Woosley &: Weaver 1982; Glatzel, El Eid, & Fricke 1985; Woosley 
1986). Alternatively, black holes with masses ^ 100 Mq may originate from an earlier population 
of stars with effectively zero metallicity (i.e.. Population III stars) and masses ^ 250 Mq (Woosley 
8z Weaver 1982; Bond, Arnett, Sz Carr 1984; Fryer, Woosley, & Heger 2001). These stars may not 
experience substantial mass loss during their lifetimes, and a massive black hole can form inside 
the star before explosive burning reverses its collapse (Bond et al. 1984). 



4.2. Accretion Disk Models 

While single component MOD and power-law continuum models have been used in the analysis 
of other IXOs (e.g., Makishima et al. 2000; Strickland et al. 2001), such an analysis may be too 
simple, leading to incorrect physical insights. The reader (and the authors) may therefore have 
misgivings about the following discussion, which is based on the results of the spectral analysis 
presented in § 3.2. 

The disk bolometric luminosity, -Lboh the inner disk radius, Rm, and the primary mass, Mx, 
are given in Table 3 for each of the spectrally modeled sources, assuming the results of the MOD 
model and that the disk is seen at = 60° [Lboi oc (cos^)^^ and oc M^ oc (cos^)~^/^, where 9 is 
the inclination of the disk (face on corresponds to ^? = 0°); see e.g., Makishima et al. (2000)]. The 
disk bolometric luminosities and color temperatures for each of the spectrally modeled sources are 
shown in Fig. 8. The lines of constant mass and Eddington ratio were added using Equations (9) 
and (11) of Makishima et al. (2000), assuming ^ = 0.412, k = 1.7, and a = 1. As can be seen in this 
figure, the majority of the sources have color temperatures similar to those observed in Galactic 
black hole candidates (i.e., /cTin ^ 1.5 keV; Tanaka 2000), primary masses of a few to a few tens of 
Mq, and luminosities which do not exceed the Eddington limit (i.e., they are consistent with being 
sources similar to the black hole candidates seen in our own Galaxy). However, Fig. 8 shows three 
sources that differ from these trends, as discussed below. 

The disk luminosities and color temperatures of the spectrally modeled sources are compared 
with the predictions of the slim disk model in Fig. 9 (cf. Fig. 1 in Watarai et al. 2001). Most of 
the sources have mass accretion rates [m = M/(LEdd/c^)] of m ^ 20, and values of Mx similar to 
those inferred from the MCD model. 
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4.2.1. J024237.9- 000118 and JO2424O.4- 000053 

The MCD model gives color temperatures of fcTin ~ 0.3 keV for two of the spectrally modeled 
sources (J024237.9-000118 and J024240.4-000053) which, together with the observed bolometric 
luminosity, yields an inner disk radius of Rm ^ 10^ km and an estimated primary mass of ^ 
100 M0 for each source (Table 3; assuming the disk to be inclined at 60°). These values are much 
larger than those derived from the spectral modeling of other IXOs (e.g., Makishima et al. 2000; 
but sec Colbert & Mushotzky 1999 for a possible exception) and Galactic black hole candidates 
(e.g., Ebisawa et al. 1994; Dotani et al. 1997; Kubota et al. 1998; Zycki, Done, & Smith 1999). 
The disk bolometric luminosities are of order 0.1 LEdd for both sources (Fig. 8), which is lower than 
that derived for other IXOs (see our Fig. 8 and Fig. 3 in Makishima et al. 2000), but consistent 
with estimated values of Lboi ^ O.lLgdd for Galactic black hole candidates in the low state (e.g., 
Nowak 1995). 

4.2.2. J024238.9- 000055 

The measured color temperature and bolometric luminosity of the most luminous source in 
NGC 1068, J024238. 9— 000055, gives an estimated primary mass of Mx ~ 2 Mq (assuming the disk 
to be inclined at 60°), which is much lower than that derived assuming the source to be radiating 
at the Eddington limit. Such low values of have been found in several other IXOs and the 
Galactic superluminal source GRS 1915+105 (e.g., Makishima et al. 2000). Mizuno et al. (1999) 
and Makishima et al. (2000) have argued that some IXOs may contain Kerr (i.e., spinning) black 
holes, so allowing a higher value of (a factor of 6 higher for a given value of and a maximally 
rotating black hole). This would not work in the case of J024238. 9— 000055, since the mass derived 
assuming the source is radiating at the Eddington luminosity exceeds the value given in Table 3 
and Fig. 8 by a factor 2± 100. Thus, we conclude that the MCD interpretation requires the source 
to be radiating above the Eddington limit, where the assumption of accretion via a geometrically 
thin disk is likely to be invalid. The slim disk model (Fig. 9) also requires J024238. 9— 000055 to be 
undergoing supercritical accretion, with m ^ 10^, which is an implausibly high accretion rate. 

The alternative description of a power-law spectrum gives a photon index of F = 0.85^[]:J^ for 
J024238. 9— 000055, which is much smaller than that seen in other IXOs and Galactic black hole 
candidates (e.g., Tanaka 2000). At low mass accretion rates, m ^ 10^'^, models of advection domi- 
nated accretion flows predict hard X-ray spectra, which are dominated by thermal bremsstrahlung 
radiation from electrons at 10^ K (e.g., Esin, McClintock, & Narayan 1997). However, there 
are two problems with such a model: (1) at energies in the Chandra band, a ~ 100 keV thermal 
bremsstrahlung continuum has a slope of F ~ 1.3—1.4, which is significantly larger than that ob- 
served in J024238. 9— 000055; and (2) such a low mass accretion rate would require a supermassive 
(^ 10^ -^0) black hole to produce the observed luminosity, which is unlikely since dynamical fric- 
tion would cause such a black hole located in the galactic disk to sink towards the center of the 
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galaxy in less than a Hubble time (e.g., Tremaine, Ostriker, &; Spitzer 1975). Strong absorption and 
hard power-law continua (0.8 ^ F ^ 1.5) characterize the 2-10 keV spectra of highly magnetized 
neutron star binary systems (e.g., White, Nagase, & Parmar 1995), and the local accretion rate 
onto the neutron star can exceed the Eddington limit by a factor of ~ 100, since the radiation 
escapes laterally out of the accretion column (Klein et al. 1996; Miller 1996). However, the bright- 
est objects of this class (e.g., the sources A0535-668, LMC X-4, and SMC X-1 in the Magellanic 
Clouds) have 2-10 keV luminosities which peak below ~ 10'^*^ erg s~^ (Liu, van Paradijs, & van den 
Heuvel 2000), an order of magnitude lower than that observed from J024238. 9— 000055. Another 
possibility is that the X-ray emission from J024238. 9— 000055 may be the result of inverse Compton 
scattering of synchroton photons in a jet (e.g., Markoff, Falcke, &: Fender 2001), an idea that could 
be tested by future coordinated radio and X-ray observations of this source. 



5. COMPARISON WITH OTHER GALAXIES 

It is of interest to compare the X-ray properties of the IXOs in galaxies with different modes 
of active star formation, e.g., circumnuclear starbursts, starbursts induced by galaxy-galaxy inter- 
actions, and star formation in "grand design" spirals. We therefore compare the luminosities of the 
X-ray point sources in NGC 1068 with those in the Antennae (Zezas et al. 2002), M82 (Matsumoto 
et al. 2001), M51 (Tcrashima & Wilson 2002a), the Circinus galaxy (Smith k Wilson 2001), NGC 
3256 (Lira et al. 2002), and MlOl (Pence et al. 2001). These galaxies have all been imaged on 
the ACIS-S3 chip, except for M82, which has been imaged with both the High Resolution Camera 
and the ACIS-I array. Unless stated otherwise, the Chandra count rates were converted to unab- 
sorbed 0.4-8 keV fluxes using PIMMS version 3. 2d (Mukai 1993) and assuming a 5 keV thermal 
bremsstrahlung continuum spectrum with absorption by the Galactic column density only. The 
values of the Galactic column density towards each of the galaxies in our sample were taken from 
Dickey & Lockman (1990). This model gives fluxes similar to those used to estimate the luminosi- 
ties in Fig. 4a for most of the spectrally unmodeled sources in NGC 1068 (cf. § 3.2). We have 
taken the average of the counts from individual sources in cases where individual galaxies have been 
observed on more than one occasion with Chandra. Following Zezas et al. (2002), a column density 
of A^H = 2 X 10^^ cm~^ was assumed for sources in the Antennae that are detected in the medium 
(2-4 keV) band, but not detected in the soft (0.3-2 keV) band (see their Table 2). Matsumoto 
&; Tsuru (1999) found strong absorption (Ah ~ 10^^ cm~^) in the residual ASCA spectrum of 
M82, obtained by subtracting the spectrum of the lowest intensity state from that of the highest 
intensity state, which may indicate the presence of strong absorption in the spectrum of the bright- 
est source in M82 (X41.4 + 60 or CXOU J095550.2-I-694047). A bright compact source in M51 
(CX0M51 J 132950. 7-1-471 155) also shows strong absorption, corresponding to Ah ~ 5 x 10^^ cm~^ 
(Terashima Sz Wilson 2002b). Therefore, intrinsic column densities of Ah = 10^^ and 5 x 10^^ cm~^ 
were assumed for these two sources, respectively. An analysis of the remaining 83 compact sources 
which project onto the disk of NGC 5194 (the main galaxy in M51) found objects that are de- 
tected in the hard (2-8 keV) band, but not detected in the soft (0.5-2 keV) band (Terashima & 
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Wilson 2002a). The hard band counts for these sources were converted to unabsorbed fluxes in 
the 0.4-8 keV band, assuming a 5 keV thermal bremsstrahlung continuum absorbed by a column 
density of A'^h = 2 x 10^"*^ cm^^ which, as noted above, is the model assumed for the sources in the 
Antennae that are not detected below 2 keV. Smith & Wilson (2001) report absorption in excess 
of the Galactic column for several of the bright compact sources in the Circinus galaxy, and so 
the 0.1-10 keV count rates were converted to fluxes in the 0.4-8 keV band assuming the results of 
their spectral modeling. The column densities from Table 3 in Lira et al. (2002) were adopted for 
the compact sources detected in NGC 3256 (we excluded the southern nucleus from our analysis 
since a low-luminosity AGN is likely to contribute to its hard X-ray emission — see Lira et al. for 
detailed discussion). 

The resulting 0.4-8 keV fluxes were used to calculate source luminosities, Lx, assuming dis- 
tances of 19.6 Mpc to the Antennae (ifo = 75 km s"^ Mpc'^; Sandage & Bedke 1994), 3.9 Mpc to 
M82 (Sakai & Madore 1999), 8.4 Mpc to M51 (Feldmeier et al. 1997), 4.2 Mpc to Circinus (Freeman 
et al. 1977), 37.6 Mpc to NGC 3256 {H^ = 75 km s'^ Mpc"^; Li'pari et al. 2000), and 7.2 Mpc to 
MlOl (Stetson et al. 1998). The total number and luminosity of sources with Lx > 2.1 x lO''^ erg s~^ 
(i.e., sources with luminosities exceeding the Eddington limit for a 1.4 Mq neutron star), observed 
42. 5-122. 5/im luminosity, star formation rate, and number of luminous X-ray sources divided by 
the star formation rate are given in Table 4 for each of the galaxies in our sample. The observed 
42.5-122.5)Um luminosity is calculated from 

Lfir = AttD^ 1.26 X 10-11 (2.58 /eo + /loo) erg s'^ (1) 

{Cataloged Galaxies and Quasars Observed in the IRAS Survey, 1985, Appendix B), where /go and 
/lOo are the IRAS flux densities in Jy at 60/xm and 100/xm, respectively, and D is the distance to 
the source in cm. We have used the flux densities given in the IRAS Faint Source Catalog (Moshir 
et al. 1990), except for M82, NGC 5194, and MlOl, for which the flux densities were taken from 
Rice et al. (1988), and Circinus, for which the flux density was taken from the IRAS Point Source 
Catalog 1988. Thermal radiation from dust heated by massive O and B type stars dominates the 
far-infrared emission between 42.5 and 122.5/Ltm in starburst and spiral galaxies (Devereux Sz Young 
1990), and so Lgi- should provide an excellent measure of the formation rate of massive stars in our 
sample of galaxies. The star formation rate (column 5 of Table 4) is calculated using Equation (26) 
in Condon (1992), and refers to stars more massive than ~ 5Mq. X-ray binaries with luminosities 
greater than 2.1 x 10^^ erg are likely to have evolved from stars with initial masses greater 
than ~ 8 Mq , so the formation rate of these stars should be similar (Grimm, Gilfanov, & Sunyaev 
2003). The Chandra observation of NGC 3256 was relatively short and so a large number of sources 
may have gone undetected in this galaxy. Further, NGC 3256 is much more distant than the other 
galaxies, so only the more luminous X-ray sources in this galaxy will have been detected. For this 
galaxy, we give a lower limit to the number of sources with luminosities exceeding 2.1 x 10^^ erg s~^. 
In fact, we note that all of the sources in NGC 3256 have unabsorbed 0.4-8 keV luminosities greater 
than 10^^ erg s"^, and so this galaxy has a larger number of IXOs than the Antennae. The galactic 
disk of MlOl at the 25.0 B-magnitude (arc sec)"^ isophote level subtends a much larger solid angle 
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than does the S3 chip, and so the number of observed sources with luminosities exceeding 2.1 x 10 

erg represents a lower limit to the actual number. We have ignored the possibility that some 
of the luminous X-ray sources may be low mass X-ray binaries; such sources would be unrelated to 
recent star formation and their presence would tend to decrease N. This effect is significant in the 
galaxies MlOl, Circinus, and NGC 5194, which have stellar masses (see e.g., Grimm et al. 2003 and 
references therein; Allen et al. 1978; Freeman et al. 1977) comparable to that of the Milky Way, 
and massive star formation rates a few times higher than that of the Milky Way (e.g., Grimm et al. 
2003). We estimate that the fraction of sources in low mass systems with luminosities exceeding 
2.1 X 10^^ erg is between ~ 25 and ~ 50% for these galaxies, assuming that the formation rate 
of such systems scales with the stellar mass of each galaxy. 

The right-hand column of Table 4 shows that, for all galaxies, except M82 and possibly NGC 
3256, the value of [N{L^ > 2.1 x 10^^ ergs-^)/S'Fi?(M0 yr"^)] = 10 to within a factor of two. 
There are considerable observational uncertainties in the measurement of both N{Lx > 2.1 x 
10^^ erg) and the SFR. The value of N is sensitive to absorption within the galaxy, and the low 
rate for M82 may be, in part, a result of absorption in this highly inclined galaxy. The low N/ SFR 
ratio in NGC 3256 is obviously related to its much greater distance than the other galaxies, with 
many X-ray sources having Lx > 2.1 x 10'^^ erg s^^ being undetected. At the moment, it seems 
premature to make any statements about variations in the ratio of the number of luminous X-ray 
sources to the star formation rate. 



6. CONCLUSIONS 

The main results of our analysis of the compact X-ray source population in the Seyfert 2 
galaxy NGC 1068 are as follows: 1) We detect a total of 84 compact sources on the S3 chip, 

of which 66 arc located within the 25.0 B-magnitude (arc sec)~^ isophotc of the galactic disk in 
NGC 1068. 2) After correcting for absorption, five of these sources have 0.4-8 keV luminosities 
greater than 10^^ erg s""*^, assuming that their emission is isotropic and that they are associated 
with NGC 1068, well in excess of the Eddington limit for a 1.4 neutron star, and so we refer 
to these sources as IXOs. If these five sources are X-ray binaries accreting with luminosities that 
are both sub-Eddington and isotropic, then the implied source masses are ^ 7-100 Mq. 3) The 
brightest compact source in NGC 1068 has a spectrum which is much harder than that found 
in Galactic black hole candidates and other IXOs, which may be the result of inverse Compton 
scattering of synchrotron photons in a jet, or other physical processes, such as complex absorption. 
4) The brightest source also shows large amplitude variability, with the flux possibly decreasing 
by a factor of two in ~ 2 ks during our Chandra observation, and by a factor of five between our 
observation and the grating observations taken just over nine months later. 5) The ratio of the 
number of sources with luminosities greater than 2.1 x 10^^ erg s~^ in the 0.4-8 keV band to the 
rate of massive (> 5Mq) star formation is the same, to within a factor of two, for NGC 1068, 
the Antennae, NGC 5194 (the main galaxy in M51), and the Circinus galaxy. This suggests that 
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the production of X-ray binaries (and SNRs) per massive star is approximately the same value for 
galaxies undergoing different types of current star formation activity. 
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Fig. 1. — The Chandra image of NGC 1068 on the scale of the galactic disk in the 0.4 to 5.0 keV 
band. The shading is proportional to the logarithm of the intensity, and ranges from (black) to 
800 (white) counts pixel^^. The linear feature running from P.A. = 78° through the nucleus to 
P.A. = 258° is a readout trail, and is not real. 

Fig. 2. — An optical image of NGC 1068 from the Digitized Sky Survey blue plate showing the 
orientation of the S3 chip (white square) with respect to the galactic disk of NGC 1068 at the 25.0 
B-magnitude (arc sec)~^ isophote level (dashed line). Also shown in this image are the locations 
of the 84 sources (crosses) discussed in § 3.1 and listed in Table 1. Coordinates are for J2000.0. 

Fig. 3. — X-ray spectra of the two brightest compact X-ray sources in NGC 1068: J024237. 9— 000118 
(filled circles) and J024238. 9— 000055 (crosses). The data are shown together with the folded model 
comprising the best-fit absorbed power-law continuum. 

Fig. 4. — Luminosity dependence of the hardness ratio (7(1.5-5.0 kcV)/C(0.4-1.5 keV) for the 
X-ray sources in NGC 1068. Horizontal dashed lines correspond to photon indices of F = — 1, 0, 1, 
2, 3, 4, and 5 when an absorbed power-law model is assumed. Filled squares correspond to sources 
for which luminosities have been estimated from the results of a spectral analysis (§ 3.2). Filled 
circles depict sources whose luminosities were estimated assuming the spectrum is an absorbed 
power-law continuum (§ 3.2). Open circles and open squares represent sources detected in only 
the hard and soft band, respectively. Three sources were not detected in the whole band (0.4- 
5.0 keV), were detected in the hard band, and were not detected in the soft band (due to the 
presence of surrounding extended soft emission). These sources are represented by open circles 
with upward and rightward pointing arrows, (a) The Galactic column density is assumed for 
sources too weak for spectral analysis, (b) The Galactic column density plus an intrinsic column 
density of A^h = 1-28 x 10^^ cm~^ is assumed for sources too weak for spectral analysis. 

Fig. 5. — The distribution of luminosities of the compact X-ray sources, assuming they are asso- 
ciated with NGC 1068. Luminosities were estimated either from the results of a spectral analysis 
or assuming the spectrum is an absorbed power-law continuum (§ 3.2). In the latter case, the 
absorbing column was taken to be either the Galactic column density towards NGC 1068 (solid 
line) or the Galactic plus an intrinsic column density of Ah = 1.28 x 10^^ cm^^ (dotted line). 
The distribution of luminosities of unrelated foreground and background sources (dashed line) was 
estimated from the spectra of sources that we detect in the field of 3C 47, assuming a power-law 
continuum model and Galactic absorption. 

Fig. 6. — The light curve for J024238. 9— 000055 binned so that the duration of each data point 
is 2000 s. The crosses represent the X-ray count rate within the source region, and include both 

source and background photons. Error bars correspond to iblo". The open triangles correspond to 
the X-ray count rate within the background region, but multiplied by the ratio of the source to 
background areas. The errors on the background rate are comparable to the size of the triangles. 



Fig. 7.— Optical continuum image of NGC 1068 at 6100 A (from Pogge k De Robertis 1993). The 
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shading is proportional to the logarithm of the intensity. Crosses mark the locations of sources 
with 0.4-8 keV luminosities (after correcting for Galactic absorption only; § 3.2) greater than 
-|^q39 gj,g g-i^ Asterisks identify sources with 0.4-8 keV luminosities (after correcting for Galactic 
absorption only) in the range 5 x 10^^ erg to 10^^ erg s~^. Open circles indicate 4 additional 
sources with luminosities (after correcting for both the Galactic column density and an intrinsic 
column density of iVn = 1.28 x 10^^ cm~^; § 3.2) greater than 5 x 10^^ erg s"^ but less than 
10^^ erg s~^. 

Fig. 8. — The disk bolometric luminosity, i^boh is plotted versus color temperature, kTi^, for each 
of the spectrally modeled sources in NGC 1068. The straight lines represent the predictions of the 

multi-color disk blackbody model of Makishima ct al. (2000). Dash-dotted lines represent the loci 
of objects with masses equal to 0.3, 3, 30, and 300 Mq. Dashed lines indicate the loci of objects 
with luminosities equal to 77 = 0.01, 0.1, 1, 10, and 100 times the Eddington luminosity, according 
to this model. 

Fig. 9. — The disk bolometric luminosity, Lboli is plotted versus color temperature, A; Tin, for each 
of the spectrally modeled sources in NGC 1068. Dash-dotted lines represent the loci of black 
holes with masses equal to 5, 10, 32, and 100 M©. Dashed lines indicate the loci of objects with 
normalized mass accretion rates equal to m = 3, 10, 32, 100, and 320, according to the predictions 
of the slim disk model of Watarai (2002). 
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50.0 ± 9.3 


0.5 ± 


0. 


,1 


38.6 


J024240.7- 


-000100 


no 
Uz 


/to 

4Z 


/in 

4U, 


vn 
. ( U 


nn 
— UU 


nn 
UU 


oy.y 


84 ± 12 


25.8 ± 9.0 


58.2 ± 9.1 


2.2 ± 


0. 


,9 


38.7 


J024240.9- 


-000144 


no 
Uz 


/lO 

4z 


A n 
4U, 


.80 


nn 
— UU 


ni 
Ui 


4o.D 


34.3 ± 7.2 


7.6 ± 4.3 


26.7 ± 6.3 


3.5 ± 


2. 


,1 


38.3 


J024240.9- 


-000027 


no 
uz 


/lO 
4Z 


/in 

4U, 


.yo 


nn 
— UU 


nn 

UU 


07 
Z ( .z 


152 ± 18 


91 ± 16 


61.8 ± 9.5 


0.7 ± 


0. 


,2 


38.8 


J024241.0- 


-000125 


no 
uz 


/lO 
4Z 


/in 

4U, 


.yo 


nn 

— UU 


ni 

Ul 


Of; 
zo.o 


38.7 ± 8.1 


24.0 ± 7.0 


14.7 ± 5.0 


0.6 ± 


0. 


,3 


38.2 


J024241.0- 


-000051 


no 
uz 


/lO 
4Z 


/in 

4U, 


Q7 

.y < 


nn 

— UU 


nn 

UU 


c;i n 

01. U 


64 ± 24 


68 ± 23 


< 7.2* 


< 0. 


.1 




38.2 


J024241.0- 


-000038 


02 


42 


40, 


.98 


—00 


00 


38.0 


80 ± 18 


56 ± 17 


24.4 ± 7.1 


0.4 ± 


0. 


,2 


38.3 


J024241.1- 


-000039 


02 


42 


41, 


.14 


-00 


00 


39.4 


136 ± 25 


99 ± 24 


36.3 ± 8.9 


0.4 ± 


0. 


,1 


38.6 


J024241.2- 


-000040 


02 


42 


41, 


.19 


-00 


00 


39.8 


126 ± 24 


89 ± 23 


37.0 ± 8.7 


0.4 ± 


0. 


,1 


38.5 


J024241.4- 


-000057 


02 


42 


41, 


.37 


-00 


00 


57.1 


9.0 ± 6.0 


< 4.7f 


12.9 ± 4.9 


> 2. 


,7 




37.7 


J024241.4- 


-000215 


02 


42 


41, 


.38 


-00 


02 


15.2 


31.3 ± 6.8 


12.6 ± 4.7 


18.6 ± 5.5 


1.5 ± 


0. 


,7 


38.1 


J024241.5- 


-000253 


02 


42 


41, 


.47 


-00 


02 


53.2 


8.7 ± 4.2 


6.0 ± 3.6 


< 6.4^ 


< 1. 


.1 




37.5 


J024241.6- 


-000123 


02 


42 


41, 


.58 


-00 


01 


22.6 


15.1 ± 5.8 


14.7 ± 5.6 


< 3.4^ 


< 0. 


.2 




37.6 


J024241.6- 


-000110 


02 


42 


41, 


.61 


-00 


01 


09.6 


12.6 ± 5.3 


6.6 ± 4.5 


6.0 ± 3.6 


0.9 ± 


0. 


,8 


37.7 


J024241. 7+000002 


02 


42 


41, 


.68 


+00 


00 


01.7 


29.6 ± 7.8 


24.3 ± 7.2 


5.3 ± 4.0 


0.2 ± 


0. 


,2 


37.9 


J024241.9- 


-000028 


02 


42 


41, 


.88 


-00 


00 


28.5 


39 ± 16 


< 44f 


12.0 ± 5.2 


> 0. 


,3 




38.0 


J024241.9- 


-000020 


02 


42 


41, 


.90 


-00 


00 


20.4 


24 ± 12 


< 29^ 


9.0 ± 4.4 


> 0. 


,3 




37.8 



Table 1 — Continued 



Name 


Combined'^ 


Soft^ 


Hard^ 




log Lx*" 


cxou 


R. A. (2000) Decl.(2000) Counts 


Counts 


Counts 


HR<i 


(erg s"^) 



J024242, 


.0-000136 


02 


42 


41 


.98 


J024242, 


.2-000047 


02 


42 


42, 


.16 


J024242, 


.5-000052 


02 


42 


42, 


.52 


J024242, 


.5-000048 


02 


42 


42, 


.54 


J024242, 


.6+000023 


02 


42 


42, 


.56 


J024242, 


.8-000245 


02 


42 


42, 


.78 


J024242, 


.8-000100 


02 


42 


42, 


.80 


J024243, 


.3-000141 


02 


42 


43, 


.34 


J024244, 


.0-000035 


02 


42 


43, 


.99 


J024244, 


.4+000435 


02 


42 


44, 


.40 


J024244, 


.6+000256 


02 


42 


44, 


.56 


J024245, 


.1-000010 


02 


42 


45, 


.12 


J024245, 


.4+000256 


02 


42 


45, 


.36 


J024245, 


.4-000219 


02 


42 


45, 


.43 


J024245, 


.5-000030 


02 


42 


45, 


.54 


J024245, 


.8+000006 


02 


42 


45, 


.79 


J024247, 


.4+000028 


02 


42 


47, 


.40 


J024247, 


.4+000351 


02 


42 


47, 


.41 


J024247, 


.5+000501 


02 


42 


47, 


.54 


J024248, 


.3+000329 


02 


42 


48, 


.32 


J024248, 


.6+000246 


02 


42 


48, 


.59 


J024250, 


.6+000546 


02 


42 


50, 


.58 


J024250, 


.9+000431 


02 


42 


50, 


.95 


J024251, 


.0-000032 


02 


42 


50, 


.99 


J024252, 


.0+000343 


02 


42 


51, 


.95 


J024252, 


.8-000055 


02 


42 


52, 


.76 


J024253, 


.1-000216 


02 


42 


53, 


.07 


J024253, 


.7+000309 


02 


42 


53, 


.70 



-00 


01 


35. 


,7 


5.2 ± 3.6 


< 4.7^ 


-00 


00 


47. 


,0 


< 16^ 


< 8.5^ 


-00 


00 


51. 


,6 


16.3 ± 6.2 


15.0 ± 6.0 


-00 


00 


48. 


,1 


24.6 ± 7.1 


< 12^ 


+00 


00 


22. 


,9 


11.2 ± 5.0 


5.5 ± 4.1 


-00 


02 


44. 


,8 


57.4 ± 


8.8 


37.9 ± 7.3 


-00 


00 


59. 


,6 


6.6 ± 4.4 


< 5.3^ 


-00 


01 


40. 


,7 


117 ± 12 


64.5 ± 9.5 


-00 


00 


35, 


,3 


92 ± 11 


42.4 ± 8.5 


+00 


04 


35, 


,1 


23.9 


± 


6.5 


7.8 ± 4.3 


+00 


02 


56. 


,1 


29.0 


± 


6.5 


16.0 ± 5.1 


-00 


00 


10. 


,0 


63.6 




9.3 


47.8 ± 8.2 


+00 


02 


56. 


,3 


16.1 




5.5 


9.8 ± 4.6 


-00 


02 


19. 


,0 


54.0 


± 


8.6 


35.3 ± 7.2 


-00 


00 


30. 


,2 


33.8 


± 


7.5 


12.1 ± 5.5 


+00 


00 


06. 


.0 


19.4 


± 


5.8 


9.2 ± 4.5 


+00 


00 


27. 


,7 


171 


± 


14 


26.7 ± 6.6 


+00 


03 


50. 


,7 


36.3 




7.7 


132 ± 13 


+00 


05 


00. 


,5 


226 


± 


17 


170 ± 14 


+00 


03 


29, 


,4 


22.9 


± 


6.4 


15.8 ± 5.3 


+00 


02 


46, 


,3 


38.9 


± 


7.9 


26.1 ± 6.8 


+00 


05 


45. 


,8 


9.1 ± 4.8 


< 2.3*^ 


+00 


04 


30. 


,6 


13.2 ± 5.2 


< 8.4^ 


-00 


00 


31. 


,6 


207 ± 16 


123 ± 12 


+00 


03 


42. 


,7 


16.1 ± 5.8 


11.4 ± 4.9 


-00 


00 


55. 


,5 


9.5 ± 4.5 


8.0 ± 4.2 


-00 


02 


15. 


,5 


20.9 ± 5.8 


17.2 ± 5.4 


+00 


03 


08. 


,7 


25.7 ± 6.5 


12.0 ± 5.0 



3.9 ± 3.2 


> 0. 


,8 


37.3 


7.4 ± 4.0 


> 0. 


,9 


37.8 


< 4.7^ 


< 0. 


,3 


37.7 


19.0 ± 5.5 


> 1. 


,6 


38.1 


5.8 ± 3.6 


1.1 ± 


1.0 


37.6 


19.6 ± 5.6 


0.5 ± 


0.2 


38.2 


5.6 ± 3.6 


> 1. 


,1 


37.4 


52.8 ± 8.4 


0.8 ± 


0.2 


38.6 


49.4 ± 8.2 


1.2 ± 


0.3 


38.5 


16.0 ± 5.5 


2.0 ± 


1.3 


38.1 


13.0 ± 4.8 


0.8 ± 


0.4 


38.0 


15.9 ± 5.1 


0.3 ± 


0.1 


38.2 


6.2 ± 3.9 


0.6 ± 


0.5 


37.8 


18.7 ± 5.5 


0.5 ± 


0.2 


38.2 


21.7 ± 5.8 


1.8 ± 


0.9 


38.3 


10.2 ± 4.5 


1.1 ± 


0.7 


37.9 


38.8 ± 7.3 


1.5 ± 


0.5 


38.8 


9.7 ± 4.8 


0.07 ± 


0.04 


37.9 


55.9 ± 9.0 


0.33 ± 


0.06 


38.8 


7.1 ± 4.3 


0.4 ± 


0.3 


37.8 


12.8 ± 5.0 


0.5 ± 


0.2 


38.0 


10.9 ± 4.7 


> 4. 


,7 


37.9 


9.1 ± 4.4 


> 1. 


,1 


37.7 


84 ± 10 


0.7 ± 


0.1 


38.9 


< 9.4^ 


< 0. 


,8 


37.8 


< 4.9^ 


< 0. 


,6 


37.6 


3.7 ± 3.2 


0.2 ± 


0.2 


37.8 


13.6 ± 4.9 


1.1 ± 


0.6 


38.0 



Table 1 — Continued 



Name 






Combined'^ 


Soft^ 


Hard'^ 


log Lx" 


cxou 


R. A. (2000) 


Decl.(2000) 


Counts 


Counts 


Counts 


HRd (erg s"^) 


J024255.9-000128 


02 42 55.92 


-00 01 27.6 


37.2 ± 7.3 


29.0 ± 6.5 


8.2 ± 4.3 


0.3 ± 0.2 38.0 



'^O. 4-5.0 keV background subtracted counts. The errors here, and in columns 5 and 6, are 
^[cr^ + (T^ X (At/v4b)^] (e.g., Bevington & Robinson 1992, p. 41), where (Tt and (Jb are the respective la 
errors on the total and background counts, derived from equation (7) in Gehrels (1986), and and Ab are 
the respective areas of the source and background regions. 

''0.4-1.5 kcV background subtracted counts. 

^1.5-5.0 keV background subtracted counts. 

•^Observed C(1.5-5.0 keV)/(7 (0.4-1.5 keV) hardness ratio. 

^0.4-5.0 keV luminosity corrected for Galactic absorption only, assuming that the source is in NGC 1068 
(distance = 14.4 Mpc) and the emission is isotropic. 

^Upper limits are calculated for a confidence level of 90%, using either a) equation (9) in Kraft, Burrows, & 

Nousek (1991) for instances when the total number of counts in the detection region is < 20, or b) equation (3) 
in Gehrels (1986), when the total number of counts is > 20. The mean number of background counts was then 
subtracted from these calculated upper limits. 



Table 2. Spectral Models of the Sources with > 50 Counts (after Background Subtraction) and 

S/N > 7. 



Source Name 

cxou 


fl022 cm-2) 




fkeVl 
ytsii V J 


Norm'l 






(d.o.f.) 


J024228.5-000202 
J024228.5-000202 


< 0.099 

< 0.094 


1 36+°-^^ 




(1.36l°;4g7) X 10-6 
(7lf ) X 10-5 


(i.06l°:^5) ^ 10-14 

X 10-15 


{2.lt\l) X 10'' 
(2.412-1) X 1038 


9.5 (3) 
11.1 (3) 


J024233.2-000105 
J024233.2-000105 


< 0.16 

< 0.10 


-. 99+0.59 
-'-■^"'-0.36 




(1.20iS:i) X 10-6 
(4lf ) X 10-5 


(i.io;°ij) X 10-14 

(9.92l°:42) X 10-15 


(2.8;?-^) X 1038 
(2.5l{:°) X 1038 


2.5 (3) 
3.2 (3) 


J024237. 7-000030 
J024237. 7-000030 


< 0.36 

< 0.17 




5-oiro 


(10+3^) X 10-^ 
(3lf 6) X 10-6 


(i.5iJ:^) X 10-14 

(1.410:3) X 10-14 


(3.8^2-9) X 1038 

(3.61°;^) X 1038 


1.4 (2) 

1.5 (2) 


J024237.9-000118 
J024237.9-000118 


n 44+0-20 

0-11+^09 

— u.uy 




0.309tg°92 

— u.uu / 


(4.4+f ^) X 10-5 

— u.o 


(2.78+°-56) X 10-14 
(2.59l°:2[) X 10-14 


(8+57) ^ 1039 
(1.1+1°) ^ 1039 


30.1 (20) 
33.7 (20) 


J024238.9-000055 
J024238.9-000055 


0.511°;}^ 
0.472^^61 


0.85i°:}? 




(4.oit°:^i) X 10-5 

X 10-5 


(5.1i;S |3) X 10-13 

(4.98l0;29) X 10-13 


(1.4631°:°^^) X 1040 
(1.408l°-°47) X 1040 


109.6 (90) 
106.8 (90) 


J024239.4-000035 
J024239. 4-000035 


0.2ll°f 

n n/s+0.30 
0-06lo.06 


■■■•^-0.8 


1 1+1.0 
1-1-0.4 


(3.3ij:^) X 10-6 

(4_4 ) X 10 


(1.24t0:63) X 10-14 

(1.00_Q 34) X 10 


(4.211:3) X 1038 
(2.8_Q g) X 10^° 


3.1 (8) 
3.0 (8) 


J024239.5-000028 
J024239.5-000028 


U.UOD_Q056 

< 0.061 


9 1 9+0.41 


0.58+°;!^ 


(5.7l4;5) X 10-6 
{Stf) X 10-3 


(2.ii;r7[;) X 10-14 

(1.38l|i;i) X 10-14 


(6.4+5-9) ^ 1038 

(3.76;°i6) X 1038 


8.0 (8) 
10.1 (8) 


J024239.7-000101 
J024239.7-000101 


r, -■ 97+0.076 
^•-•-^'-0.080 

< 0.044 


i-93lg:i 


0.981°:?° 


(1.96l°:69) X 10-5 
{2>At\l) X 10-3 


(7.6li-^) X 10-14 
(6.0lj;i) X 10-14 


(2.42+°-4i) X 1039 

(1.56l°;26) X 1039 


31.2 (32) 
33.0 (32) 


J024240.4-000053 
J024240.4-000053 


n 70+0.30 
'J- '^-0.21 
S5+0-20 


5.8l}i 


r, 970+0.076 
U.Z(O_Q060 


{Q.^tll) X 10-5 
77+4-21 

'-0.61 


(1.53l0i6) X 10-14 
(1.47;°i2) X 10-14 


(2+32) ^ ^q40 
(1.412.2) X 1039 


22.9 (18) 
23.8 (18) 



Table 2 — Continued 



cxou 


(10^2 cm-2) 




(keV) 


Norm'l 


(erg cm-2 s-^) 


(erg s-i) 


fd f ) 


J024240.5-000037 
J024240.5-000037 


< 0.10 


9 (-7+0.53 


0.48l°:i^ 


(6.0l|^) X 10-6 


(1.59l°-i) X 10-^' 
(1.161°:!^) X 10-14 


(6lf ) X 1038 
(3.211:6) X 1038 


7.9 (14) 
9.2 (14) 


J024240.7-000055 
J024240. 7-000055 


< 0.14 

< 0.067 


1 47+0.41 

—0.28 




(4.7lf J) X 10-6 
(3.4l|9) X 10-4 


(3.2710-50) X 10-14 
(2.66lo;5^0) X 10-14 


(8.451J-15) X 1038 
(6.8lJ:I) X 1038 


32.2 (20) 
35.1 (20) 


J024240.7-000100 
J024240.7-000100 


0.901J-38 
0.601H? 


1 ^q+1.22 


3+00 


(5+f ) X 10-6 
(51^0) X 10-5 


(2.6610-81) X 10-14 
(2.5210;4^) X 10-14 


(lOlf) X 1038 

(7.711:8) X 1038 


13.9 (14) 
14.3 (14) 


J024240.9-000027 
J024240.9-000027 


< 0.092 

< 0.043 


1-661^:28 


±.uo_Q 28 


(5.0+f 6) X 10-6 
(l.Ol^-O) X 10-3 


(2.87l0-5g^) X 10-14 
(2.10111:62) X 10-14 


(7.5 ± 1.3) X 1038 
(5.5lJ:5) X 1038 


14.7 (16) 
19.7 (16) 


J024242.8-000245 
J024242.8-000245 


< 0.16 

< 0.064 


i.59i°:l? 


1 1+1-3 


(1.5lii) X 10-6 
(2lf ) X 10-4 


(9.211^) X 10-15 
(7.0l4|) X 10-15 


(2.39 ± 0.76) X 1038 

(1.811:?) X 1038 


1.2 (3) 

5.3 (3) 


J024243.3-000141 
J024243.3-000141 


0.09l°:J3 
< 0.14 


l-49lSi8 


1 4+1-3 


(3.2t3-i) X 10-6 
{2tl^) X 10-4 


(2.01J-1) X 10-14 

(1.611:^) X 10-14 


(5.6l?-2) X 1038 

(4.01^:9) X 1038 


2.3 (4) 
2.6 (4) 


J 024244. U-UU0U35 
J024244.0-000035 


r\ /ioH~0.33 

0-42_o.i9 

ri r)f)+0.25 
^•^^-0.15 


1 nA+0.72 


-, 90+0.62 
-'--^'*-0.36 


(5.5^2 6/ X 10 
(4113) X 10-4 


(1.74Iq g-^ j X 10 
(1.53l°;3?) X 10-14 


o+6.4\ , , 1 n38 

(d.8_i 3) X 10 
4.7l°:?) X 1038 


>'7\ 

3.0 (7) 
2.4 (7) 


J024245. 1-000010 
J024245. 1-000010 


< 0.15 

< 0.089 


2 0+1-'' 


0.53l°:48 


(1.8ll:I) X 10-6 
(4lf ) X 10-3 


(8.011;^) X 10-15 
(4.714-5) X 10-15 


(2.1+3-8 ^ 1038 
(1.311:1) X 1038 


1.3 (3) 
4.2 (3) 


J024245.4-000219 
J024245.4-000219 


26+°-34 

n5+°-22 


/iq+0.65 
^•^^-0.93 


0-86t°i^ 


(2.814-5) X 10-6 
(6lf ) X 10-4 


(6.311-1) X 10-15 
(5.6l|i) X 10-15 


(3+P) X 1038 

(1.5810:6^) X 1038 


0.51 (2) 
0.72 (2) 



Table 2 — Continued 



Sim 1 TT'P l\l aTYi 

cxou 






(keV) 


Norm*^ 


^ X 

(erg cm-2 s-^^ 


) 


(erg s-^) 


v2 (d nf) 
X (,1-1. u.i.; 


J024247.4+000028 
J024247.4+000028 


"■"-'■^'-'—0.028 

< 0.032 


9 1 y+0.42 
' —0.36 


o.6ol°:l« 


(5 0+?-°) X 10-6 
(6lf ) X 10-3 


(1 89+n1o) X 10' 
V-^-""^— 0.49/ 

(1.35l°:2f) X 10- 


-14 
-14 


(5 5+3-2) JQ38 

(3.65l°i?) X 1038 


2 7 (8) 
9.7 (8) 


J024247.5+000501 
J024247.5+000501 


< 0.067 

< 0.024 


^•^^-0.24 




^•^^-0.14 


(6.1+J-9) X 10-6 
(l.llgi) X 10-2 


(2.641°:^,) X 10- 
(IMt'oil) X 10- 


-14 
-14 


(7.1+J-3) X 1038 

X 1038 


5.1 (12) 
21.6 (12) 


J02425 1.0-000032 
J024251. 0-000032 


4+0.078 

< 0.034 


J^-^^-0.25 


1 28+°-^^ 


(5.212-4) X 10-6 
(5.7l|3) X 10-^ 


(3.9l°-^) X 10- 
(2.92lg:?2) X 10- 


14 

-14 


{IMt'oiD X 1039 
(7.5l?:3) X 1038 


4.8 (10) 
6.3 (10) 



Note. — The errors given here are at the 90% confidence level for one interesting parameter (i.e., Ax^ = 2.706). 
^Measured intrinsic column density. 
^Photon spectral index. 

'^Inner disk temperature for multi-color disk blackbody model. 

'^Normalizations for the power-law and MCD models. For the power-law model, the units are photons cm-2 s"^ keV"^ at 1 keV. 
For the MCD model, the normalization is equal to [(rin/km)/(D/10 kpc)]2 cos0, where nn is the inner disk radius, D is the distance 
to the source, and 9 is the inclination of the disk (face on corresponds to = 0°). 

^Observed 0.4-8 keV flux. 

^0.4-8 keV luminosity corrected for both Galactic and intrinsic absorption, assuming that the source is in NGC 1068 (assumed 
distance = 14.4 Mpc). Luminosities are calculated assuming the emission to be isotropic for the power-law model, and the disk to 
be inclined at 60° for the MCD model. 



Table 3. MCD model parameters. 
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Note. — Quantities in the table scale as L^oi oc 
(cos0)~-^ and Mx oc oc (cos0)~^/^, where 6 is the 
inclination of the disk. We have assumed cos 9 = 0.5, 
corresponding to the angle averaged value, = 1.18, 
and a = 1 (§ 4.2). 

^Disk bolometric luminosity calculated using equation 
(3) in Makishima et al. (2000). 

^Inner disk radius, derived from equation (5) in Mak- 
ishima et al. (2000). 

^Mass of the compact object estimated from equation 
(8) in Makishima et al. (2000). 



Table 4. Luminosities of the Compact X-ray Sources in Nearby Galaxies. 



Galaxy 



N{L^ > 2.1 X 10^^ erg 3-^)=^ 



log (L 

X, total 

log (Lfii-^ 
(erg s~^) (erg s~^^ 



SFR'^ 
(Moyr-i) 



Af(Lx>2.1xlO^» ergs-^) " 
SFR 



NGC 1068 24 40.2 44.3 5.0 5.0 

the Antennae 28 40.6 44.0 2.4 12 

M82 6 40.7 44.0 2.2 2.8 

NGC 5194 11 40.1 43.6 0.9 12 

Circinus 9 40.4 43.4 0.6 15 

NGC 3256^ > 13 > 41.0 44.9 17.4 > 0.7 

MlOl^ >5 > 39.6 43.6 0.9 > 5.6 



^Number of sources with Lx > 2.1 x 10^^ erg s ^ in the 0.4-8 keV band. 

'^Total luminosity from sources with Lx > 2.1 x 10^^ erg in the 0.4-8 keV band. 

■^Observed 42. 5-122. 5^m luminosity (see § 5). 
•^Star formation rate (see § 5). 

^The limiting luminosity level during the Chandra observation was ~ 10^^ erg s"^ in the central region of this 

galaxy. 

^The galactic disk of MlOl at the 25.0 B-magnitude (arc sec)~^ isophote level subtends a much larger solid 
angle than does the S3 chip. 
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